Liesegang patterns: Effect of dissociation of the invading electrolyte 
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The effect of dissociation of the invading electrolyte on the formation of Liesegang bands is 
investigated. We find, using organic compounds with known dissociation constants, that the spacing 
coefficient, 1 + p, that characterizes the position of the n-th band as x n ~ (1 + p)", decreases 
with increasing dissociation constant, Kd- Theoretical arguments are developed to explain these 
experimental findings and to calculate explicitly the Kd dependence of 1 + p. 



I. INTRODUCTION 

Liesegang patterns are quasi-periodic structures pre- 
cipitated in the wake of a moving reaction front. Al- 
though they appear only in specific physico-chemical con- 
ditions, these structures are widespread in nature and can 
be found in systems ranging from biological (populations 
of bacteria) to geological (structures in agate rocks) . 
Depending on the geometry and the dimensionality of the 
system, the observed patterns are bands, rings or spheres, 
although it is also possible to generate more particular 
patterns such as spirals. 

In chemistry, such structures are produced by allow- 
ing two chemicals, called A and B, initially separated, 
to mix through diffusion. The two components are cho- 
sen so that they react and form a nonsoluble product, 
called precipitate. In a typical example, B (the inner 
electrolyte) is initially dissolved in a gel and placed in a 
test tube. At time t — 0, the reagent A (outer electrolyte) 
is poured into the tube and it starts to diffuse into the 
gel. As the result of chemical and physical mechanisms 
such as nucleation, aggregation, coagulation or floccula- 
tion, opaque or colored zones called Liesegang structures 
are formed as the diffusive front of A moves ahead (see 
Fig.l for an example). 

The precipitation zones that appear in the tube usu- 
ally exhibit a few well defined features. In particular, 
the time of the formation of the n-th band, t n , and the 
distance of the n-th band from the initial interface, x n , 
are related by the so called time law Q : 

Xn ^° \/^n (1) 

This law follows evidently from the diffusive nature of 
the reaction-diffusion front. 

A more intricate property of the positions of the bands 
is that they usually form a geometric series, x n ~ (l+p) n , 
implying that x n +x/x n , tends to a constant for large n. 
This experimental observation is usually referred to as 
the Jablczynski law or spacing law H : 



X n +1 



The quantity 1 + p will hereafter be referred to as the 
spacing coefficient. Usually p > and one speaks of reg- 
ular banding. 

A central goal of the studies of Liesegang bands is to 
understand how the spacing coefficient depends on both 
the experimentally controllable parameters (such as the 
initial concentrations of A and B) and the other less con- 
trollable material parameters such as the diffusion coeffi- 
cients. A major achievement from experimental point of 
view has been the establishment of the so called Matalon- 
Packter law ]|,|| ■ It gives the dependence of the spacing 
coefficient on the initial concentration, ao, of the outer 
electrolyte, A, in a simple form 



l+P = Qi 



02 

«o 



(3) 



where Q\ and Q2 depend on the initial concentration, &o , 
of the inner electrolyte as well as on all the other material 
parameters. 

As far as the theoretical approaches are concerned, 
several competing theories have been developed [0-^8), 
many of whic h M |l4| fare well in deriving the time and 
spacing laws [[19|. They succeed in explaining the forma- 
tion of these bands by following how the diffusive reagents 
A and B turn into immobile precipitate D 



A + B^ ...C... 



D 



(4) 



= 1 



n>l 1 
Pn ► 1 



P 



(2) 



taking into account various scenarios for the intermedi- 
ate steps denoted here by ...C... . The theories are not 
entirely equivalent when they are applied to derive the 
Matalon-Packter law and comparison with experiments 
slightly favors the so called induced sol- coagulation model 
p0| ]. It should be emphasized, however, that there is no 
clear experimental evidence at present that would dis- 
tinguish decisively among the existing theories. Accord- 
ingly, it remains to be a task to search for distinguish- 
ing features of Liesegang phenomena both in experiments 
and theories. 
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From a chemical point of view the process given by 
(||) often means oversimplification. In some cases, for 
example, the precipitate D may partly redissolve in the 
surplus of component A. This is not a negligible effect 
since the precipitation and dissolution processes, coupled 
to diffusion, may lead to a propagating band instead of 
a final static pattern . 

Another effect, that is our main concern in this paper, 
is the dissociation of the invading electrolyte A. In the 
theoretical approaches formulated up to now, it is usu- 
ally assumed that the concentration of the reacting ions 
is equal to the concentration of the electrolyte. This is 
true only when an ionic compound, a strong acid or a 
strong base is used as A. If, however, a weak acid or 
weak base is used as outer electrolyte, then the two con- 
centrations are not equal and the effect of dissociation 
must be taken into account. 

As an explicit example, imagine the formation of 
Mg(OH)2 precipitate (this process will be our main con- 
cern throughout this paper). Consider a gel soaked with 



MgCl2, and let us pour either NaOH or NH4OH onto 
its surface. In both cases the chemistry is the same: 
Mg(OH)2 precipitate is formed. There are, however, 
significant differences in the patterns formed as shown in 
Fig.l. The strong base (NaOH) results in a very dense 
band structure, whereas the weak base (NH4OH) pro- 
duces a rather loose pattern. A more precise analysis 
reveals (see below) that the latter structure has a signifi- 
cantly larger spacing coefficient. Based on this example, 
one can conclude that either the co-ions or the different 
degree of dissociation have a strong influence on the final 
pattern. In previous experiments on Liesegang phenom- 
ena, we found that the effects of co-ions (such as Na + 
or NH^) is almost negligible [Q. We shall therefore as- 
sume that the decisive role is played by the dissociation 
and that it may be possible to isolate the effect of disso- 
ciation on the formation of Lieasegang bands. Since we 
aim at studying this effect and since we count physicists 
in our audience, a few words are in order about dissoci- 
ation. 




FIG. 1. Mg(OH)2 Liesegang patterns formed from the reaction of the inner electrolyte MgCh with two different outer 
electrolytes: NH4OH for the left picture and NaOH for the right one. 



Dissociation is a common phenomenon for acids and 
bases. Strong acids and bases dissociate into ions en- 
tirely i.e. the solution contains only ions but no neutral 
molecules. For weak acids and bases the dissociation is 
not complete and the neutral molecules and ions keep a 



dynamical equilibrium in the solution. 

In order to introduce notation, let us consider a weak 
base WOH . In an aqueous solution, the base, its cations, 
W + , and anions, OH~ , are all present 



2 



WOH 



OH- 



(5) 



pK b 



log {K d /e ) 



(10) 



and the equilibrium concentrations are related by 
[W+][OH-] 



K d 



[WOH] 



(6) 



where Kd is the dissociation constant. It is convenient to 
introduce the degree of dissociation, a defined as 



[OH- 



(7) 



a a 

where ao denotes the initial concentration of the base, 
ao = \WOH] t —Q. Then eq.(0) can be written as 



a 



1 — a 



(8) 



and one can see that a can be determined once ao and 
Kd are known. 

In order to study the effect of dissociation on Liesegang 
pattern formed by Mg(OH)2 precipitate, several bases 
having different values of Kd are required. Unfortu- 
nately, there are only few inorganic ones available. Con- 
sequently, we have chosen organic compounds to provide 
the necessary OH~ ions. The organic molecules we shall 
use dissociate like the 1-buthylamin: 



C 4 H 9 - NH 2 + H 2 ^ C 4 H 9 - NH+ + OH~ 



(9) 



For our purposes, an important characteristic of an amine 
is its basicity that is conveniently expressed by the pK a 
of its conjugate acid. For the sake of simplicity, however, 
we shall characterize the dissociation ability of amines 
by the basicity constant pK b — 14 — pK a that can be 
expressed through Kd as 



where eo = lmol/dm 3 and Kd is given by eq. (||). The 
basicity constant is essentially a measure of an amine's 
ability to accept a proton from water according to (|^). 
The higher the pK b of an amine the lower is its dissoci- 
ation ability. 

The goal of this work can be now stated more precisely. 
We carry out a combined experimental and theoretical 
study of the effect of dissociation on the Liesegang pat- 
terns by investigating the pi^-dependence of the spacing 
coefficient 1 + p. 

The paper is organized as follows. In Sec. 2, we describe 
the experimental procedures and present the experimen- 
tal data. The main finding here is that increasing the 
dissociation constant results in a decreasing spacing co- 
efficient. In Sec. 3, the theoretical models developed for 
explaining the formation of Liesegang patterns are gen- 
eralized by including the dissociation effects. Then the 
spacing coefficient as a function of the dissociation con- 
stant is calculated explicitly and comparison with the 
experimental data follows. The conclusion is that all 
theoretical models that produce the Matalon-Packter law 
should be in qualitative agreement with the experiments. 



II. EXPERIMENTAL PART 
A. Materials 

We have studied the formation of Mg(OH)2 precip- 
itate, by diffusion of inorganic and organic compounds 
into chemically cross-linked polyvinylalcohol (PVA) hy- 
drogel containing the MgCl 2 inner electrolyte. 



Outer electrolyte 


pK b 


1+p 


NH 4 OH 


4.75 


1.76 ± 0.05 


2-buthylamine 


3.44 


1.13 ± 0.01 


i-propylamine 


3.40 


1.110 ± 0.005 


1 ,2-ethyldiamine 


3.29 


1.09 ± 0.01 


1-buthylamine 


3.23 


1.09 ± 0.01 


piperidine 


2.88 


1.043 ± 0.004 


NaOH 


— oo 


1.00 ± 0.01 



Table I. Outer electrolytes, their basicity constants, and the corresponding spacing coefficients obtained from the 
experiments. Note that the case of NaOH appears to be special in that x n +i/x n — > 1 for large n and the asymptotic 
form of x n should be different from x n ~ (1 +p) n . 



Table I lists the bases used for producing the ions OH~ 
as well as their basicity constants p3| . When the spa- 
tially separated reactants come into contact, the white 
precipitate Mg(OH)2 is formed as a result of the follow- 
ing reaction: 



Mg 2 



20H- -► Mg{OH)i 



(11) 



B. PVA gels as reaction-diffusion media 



The highly swollen PVA-hydrogels were prepared by 
cross-linking of primary PVA-chains with glutaric alde- 
hyde (GDA) in aqueous solution, containing dissolved 
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MgCh- Commercial PVA (Merck 821038) and solution 
of 25 mass% GDA (Merck) were used for preparation. 
The initial polymer concentration, as well as the cross- 
linking density has been kept constant. The polymer 
content of the PVA gels was in every case 3.5 mass%. 
The ratio of monomer unit of PVA (VA) to the cross- 
linking agent, GDA was: (VA)/(GDA)= 250. In or- 
der to study the precipitation and band formation in 
swollen networks, one of the reactants, the inner elec- 
trolyte MgCh was mixed with the polymer solution con- 
taining the cross-linking agent. The concentration of 
MgCl% in the gel phase was constant (0.05mol/dm 3 ). 
The gelation was induced by decreasing the pH of the sys- 
tem by HCl (Reanal, Hungary). Then the solution was 
poured into glass tubes. The gelling process usually took 
3-4 hours. Experiments were made in tubes of length 300 
mm and diameter 12 mm. The tubes were sealed and al- 
lowed to stand undisturbed at constant temperature of 
17°C. After completion of the network formation, the 
outer reactants were brought into contact with the gels, 
so that they start to diffuse into it. The diffusion took 
place vertically downward at controlled temperature. We 
maintained the constancy of boundary condition by re- 



freshing the upper solution continuously, thereby keeping 
the concentration of the outer electrolyte at a fixed value 
(7 mol/1). During all the experiment process, constant 
temperature of 17.00 ± 0.02°C was maintained and the 
tubes were kept free from mechanical disturbances such 
as vibrations. 



C. Liesegang band formation 

The penetration of organic as well as inorganic 
molecules into the gel containing the Mg 2+ ions resulted 
in precipitate bands (see Fig. 1) with a sharp interfaces of 
Mg(OH)2- In order to determine the position of bands 
we used a digital video system. A CCD camera with 
an 1/3" video chip has been connected to a PC through 
a real-time video digitizer card. We have also used a 
scanner to digitalize the experimental results. We have 
determined the coordinate of the gel surface where the 
diffusion started. Then the position of the n-th band, x ni 
was measured as the distance between the surface of the 
gel and the upper side of the n-th precipitate. 




FIG. 2. Dotormination of the spacing coefficient given by slope of the curves x„+i vs. x n for different compounds. The 
distances are measured in pixels of the digitalized pictures. The choice of the unit length is irrelevant as we are interested in 
the ratio x n+ i/x n . The straight lines are least-square fits. 



In order to determine the value of the spacing coeffi- 
cient we plotted x n +i against x n . For the sake of a better 
visibility, only three samples are shown on Fig. 2 but the 
other sets of data are of similar quality. The data were 
analyzed by linear least square method and the resulting 
spacing coefficients and related errors are summarized in 
Table I. Note that the relatively large number of bands 



allowed a rather accurate determination of the spacing 
coefficients (of course, one should keep in mind that the 
errors quoted in Table I are the statistical errors and they 
do not include possible systematic errors). 

The pKi, dependence of the spacing coefficient is shown 
in Fig. 3. One can easily see that increasing the dissocia- 
tion (pKb = 10~ Kd ) decreases the spacing coefficient. It 
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should be pointed out, however, that in contrast to the 
theories discussed below, where all parameters except K c i 
can be kept fixed, the situation is more complex in the 
experiments. All the parameters such as reaction thresh- 
olds and diffusion constants may, in principle, change 
when going from one outer electrolyte to the other. Thus 



Fig. 3 displays the coupling between the spacing coeffi- 
cient and the value of pKb under the assumption that 
one can neglect the effect of the co-ions on the material 
parameters that are relevant in the process of pattern 
formation. 




FIG. 3. Spacing coefficient 1+pasa function of the basicity pKt of the invading electrolytes. The full line is the prediction 
of the theory, eq.(22), with parameters Qi = 1 and Q2 ~ 0.003. 



III. THEORETICAL APPROACHES 

The formation of precipitation bands is a rather com- 
plex phenomenon and it is by no means certain that there 
exist a unified description of all Liesegang phenomena. 
Indeed, several theories have been developed. All of them 
follow how the diffusive reagents A and B turn into im- 
mobile precipitate D through some intermediate steps 
A + B — *■ ...C... — ► D, The common feature of the the- 
ories is that the precipitate appears as the system goes 
through some supersaturation and nucleation thresholds. 
The differences arise in the details of treating the thresh- 
olds and in the growth kinetics of the precipitate. 

Our aim here is not to review these theories in detail (a 
recent critical discussion can be found in ^0|). Instead, 
we shall recall them briefly, highlight their specific fea- 
tures and then discuss how dissociation enters into the 
problem. We shall see that the experimentally relevant 
regime of dissociation that is fast compared to diffusional 
relaxation, can be treated independently of the details of 
the theories. The final result is that, since the Matalon- 
Packter follows from all of the theories, the consequence 
of the dissociation is universal. Namely, the spacing co- 



efficient depends exponentially on the basicity constant, 
P K b . 



A. Supersaturation of ion-product theory 



The simplest theory is based on the concept of su- 



and has been developed 
In this theory, the A, B 



per saturation of ion-product 
by many researchers ] ^Io| , |2^ 
reagents turn into the precipitate D (without any inter- 
mediate steps) provided the local product of concentra- 
tions of the reactants, ab, reaches some critical value, q* . 
The nucleated particles grow, deplete A and B in their 
surroundings and, as a consequence, nucleation stops. 
When the reaction zone (where ab is maximum) moves 
far enough so that the depletion effect of the precipitate 
becomes weak, then ab = q* is reached again and nucle- 
ation can occur. This quasiperiodic process leads to the 
formation of successive bands. In the limit of very large 
precipitation- and aggregation rates, the control param- 
eters of this model are the initial densities of the elec- 
trolytes ao and bo, their diffusion coefficients D a and Db 
and the ion-product threshold q* [M. 
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B. Nucleation-and-growth theory 

The nucleation-and-growth theory introduces a single 
intermediate step in ...C... with a mechanism of band 
formation based on the supersaturation of the interme- 
diate compound C In this theory, A and B re- 
act to produce a new diffusing species C the nature of 
which is not really specified. It may be a molecule as 
well as a colloid particle. The main event is the nucle- 
ation that occurs when the local concentration of C-s 
reaches some threshold value. The nucleated particles 
(-D-s) act as aggregation seeds and the nearby C-s aggre- 
gate to the existing droplet (hence become D-s) provided 
their local concentration is larger than a given aggrega- 
tion threshold. These models are characterized by two 
thresholds, c* for nucleation and g* for droplet growth. 
As before, ao, bo, D a , D are control parameters but the 
diffusion constant D c of the C species and c* , g* appear 
as extra control-parameters [Q. The depletion mecha- 
nism around an existing precipitation band is similar to 
the one described for the ion-product theory and it leads 
to the quasiperiodic band formation. 



C. Induced sol-coagulation theory 

The nucleation-and-growth theories have been modi- 
fied to take into account effects of the concentration of 
the electrolytes on the nucleation processes. In the in- 
duced sol- coagulation theory |]ll| , |l2| , it is assumed that A 
and B react to produce a sol (C) and this sol coagulates 
if both the concentration of C exceeds a supersatura- 
tion threshold c > c* and the local concentration of the 
outer electrolyte is also above a threshold a > a*. The 
quantity a* is often referred to as the critical coagulation 
concentration threshold and is a new free parameter in 
this theory. The band formation is a consequence of the 
nucleation and growth of the precipitate combined with 
the motion of the front where a = a*. 



D. Reaction-diffusion equations 

All the above theories can be described in terms of re- 
action diffusion equations for the concentrations of the 
reagents (a, 6, c) and of the precipitate (d): 



dta = D a d 2 a — i?i(a, b, c, d) , 
d t b = D b d 2 J) - Ri(a, b, c, d) , 
d t c = D c d 2 c + i?2(a, b, c, d) , 



d t d = 



R 3 (a,b,c, d) 



(12) 
(13) 
(14) 
(15) 



where the reaction terms R a can always be chosen so 
that they describe the precipitation and aggregation pro- 
cesses which are building blocks of the theories discussed 
above. The time- and the spacing laws (0,0) follow from 



the equations corresponding to the above theories. Fur- 
thermore, using these theories, one can also derive EG] 
the experimentally observed Matalon-Packter law (Q), 
1 +p = Qx + Qil a§. This is an important law because ao 
is one of the few experimentally controllable parameter 
in Liesegang phenomena. As we shall see, the knowledge 
of the simple ao dependence of 1 + p allows us to derive 
the effects of the dissociation on the spacing coefficient. 



E. Role of dissociation 

The outer electrolyte A produces the reagent A 
through reversible dissociation A <-> A + A'. This process 
of dissociation can be characterized by a relaxation time, 
Tdis, defined as the typical time taken by a molecule to 
dissociate and thus to equilibrate the local ionic concen- 
trations. This is a microscopic time and thus it should 
be much smaller than the time of diffusional relaxation 
of density profiles. Indeed, the diffusive front of the A 
particles moves with a velocity vt ~ W D a /t that dimin- 
ishes with time. Thus the relaxation time of density per- 
turbations over a characteristic distance I (£ can be the 
width of the reaction zone, w ~ i 1 / 6 p5| , or the distance 
between consecutive bands x n +\ — x n ~ \ft^) increases 
without bound i.e. it is a macroscopic time 



t° , cr> 1/2 



(16) 



Consequently, in the long time limit, Tdis becomes neg- 
ligible and we can assume the existence of a local dis- 
sociation equilibrium. Thus we can extend relation ^) 
to the out-of-equilibrium situation, where concentrations 
depend on x and t. This implies that, denoting the den- 
sity of the reagent ion of the outer electrolyte A by a, the 
first reaction diffusion equation ( |l2| ) must be replaced by 
the following couple of equations: 

dta = D a d 2 a — ma + (17) 
dta = Dad 2 a + ma — K2a 2 — R\(a, b, c, d) . (18) 

where ki,K2 — > 00 accounting for Tdis 3* Tdiff and 
K1//C2 = Kd ensuring that the steady, homogeneous state 
satisfies the steady state condition (^). The rest of the 
equations ( p^|Jll|Jl5| ) change only by a — > a in the reaction 
terms. 

Equations (J17|,|lq) together with the rest of the 



reaction-diffusion equations (|l3|,y_4 15) can be solved nu- 



merically and the spacing coefficient can be determined 
within the framework of the various theories. We do not 
have to carry out this work, however, since the concentra- 
tion of the outer electrolyte in a usual experimental setup 
is two orders of magnitude larger than that of the inner 
electrolyte. This means that, for practical purposes, the 
precipitation processes do not influence the concentration 
profile of A. In turn, this means that the concentration 
of A is determined by that of A and, using the fact that 
the local dissociation equilibrium is established fast, we 
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return to the original problem of equations ( |L2] , |L^ , [i4]Jl5| ) 
but with the initial concentration, ao, replaced by So as 
determined from the steady state condition, eq.(@): 



(19) 



It follows then that the effect of dissociation on the 
Matalon-Packer law can be obtained by just replacing 
a by a in eq.(g): 



1 



P = 



2Q 2 



4a Q K d - K d 



Qi 



2-10P K >Q 2 /e 



y/T+T- 10P^a /e - 1 



(20) 



(21) 



where we have replaced Kd by the experimentally mea- 
sured basicity parameter pKb = — \og(Kd/eo) in the sec- 
ond equation. This expression of the spacing coefficient 
through the basicity parameter is the central result of our 
theoretical discussion. The derivation has no recourse to 
the details of the underlying theories of precipitation and 
so the result is valid for any theory that reproduces the 
Matalon-Packter law. 

In the experiments discussed in previous chapters, the 
basicity parameter is large (pKf, ^3 — 5) and since 
a 7eo, we have 10 pKb a /eo ^ 1. As a consequence, 
the Matalon-Packter law simplifies to 



l+P = Qi 



Q-2 



(22) 



thus resulting in a simple exponential dependence on the 
basicity constant, pK^- 



F. Comparison between experiments and theory 



and the surprise here is that these dependences are not 
significant or absent entirely. The agreement actually 
points to the correctness of our initial assumption that 
the Mg 2+ and OH~ ions dominate the process with the 
co-ions playing no role at all except for the co-ion of OH~ 
setting the initial concentration of the OH~ ions. 

As we can see, the above theoretical studies of dis- 
sociation and the comparison with experiments did not 
help in deciding which was the right theory for Liesegang 
phenomena. On the other hand, these studies lead to a 
simple picture of how the dissociation affects the spacing 
coefficient and, furthermore, it has been possible to ex- 
press the results in simple analytical form that should be 
useful when discussing and designing Liesegang patterns. 



IV. SUMMARY 

Our experiments show that the spacing coefficient of 
Liesegang patterns is strongly influenced by the degree of 
dissociation of the outer electrolyte. The dependence of 
l+pon the basicity constant, pK^ has been found to be 
exponential and this experimental finding has been ex- 
plained on the basis of the present theories of Liesegang 
phenomena. 
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The functional forms obtained above (eqs.|T 22) are 
universal in the sense that the pK^, dependence is explicit 
and only Q\ and Qi depend on the details of the theories. 
Unfortunately, Q\ an d Qi contain unknown parameters 
such as aggregation thresholds and so they should be con- 
sidered as fitting parameters. Nevertheless, it remains a 
question whether the experimental data could be fitted 
with reasonable values of Q\ and Qi (the only obvious 
restriction on Q\ and Qi coming from theories is that 
Qi > 1 and Q 2 > 0). 

Fig. 3 shows the experimental data together with the 
theoretical curve ( |22|) using Q\ = 1 and Q 2 = 0.003 and 
the agreement we see is very good. This is somewhat sur- 
prising. Indeed, using various outer electrolytes changes 
not only the dissociation constants but may alter the dif- 
fusion constants of the reagents (Mg 2+ , OH~) as well as 
it may change the various precipitation thresholds. The 
spacing coefficient should depend on all those quantities 
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